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As usual, let R, C be the real and complex numbers, and for each positive 
integer n let R", C" be the real and complex vector spaces of ra-tuples of real 
and complex numbers, respectively. We write (^-^(R"), GL{C^) denote the 
general linear groups of invertible real and complex- linear transformations 
on R", C", respectively, and ^L(R"), ^L(C") for the subgroups of ^^(R"), 
GL(C") consisting of linear transformations with determinant equal to 1. 
Also, 0(R"), [/(C^) are the subgroups of ^^(R"), GL(C") of orthogonal 
and unitary linear transformations on R", C", respectively, which are the lin- 
ear transformations which preserve the standard Euclidean norm, or, equiv- 
alently, the standard inner products on R", C", which is the same as saying 
that they are invertible linear transformations whose inverses are equal to 
their adjoints, and S'0(R"'), SU{W) are the special orthogonal and unitary 
groups, which are the subgroups of 0(R"), [/(C") of linear transformations 
which also have determinant equal to 1. 

Let us write iS(R"'), iS(C"') for the real vector spaces of self-adjoint linear 
transformations on R", C". We write iS+(R"'), 5+(C"') for the open con- 
vex cones in iS(R"'), 5(C") consisting of linear transformations which are 
positive-definite. Furthermore, we write A^(R"'), M.{C"') for the subsets of 
5+(R"), iS+(C"') of linear transformations with determinant equal to 1. 

Of course GL(R"), G'L(C") are open subsets of the vector spaces of all 
linear transformations on R", C", respectively, which can be defined by 
the condition that the determinant is nonzero. One can think of S'L(R"'), 
S'L(C") as hypersurfaces in GL(R"), G'L(C"), and indeed they are regular 
or smooth hypersurfaces because the equation that the determinant be equal 
to 1 is nondegenerate on these subsets, which is to say that the gradient of 
the determinant, as a function on the vector space of linear transformations 
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on R" or C", is nonzero at every invertible linear transformation. Similarly, 
A^(R"), 7Vf(C") are smooth hypersurfaces in the open sets 5_|_(R"), 5+(C"). 
There are natural actions of G'L(R"), G'L(C") on 5(R"), 5(C"), defined 

by 

(1) A^T*AT 

when A is a self-adjoint linear transformation on R" or C" and T is an 
invertible linear transformation on the same space. If T is an invertible linear 
mapping on R" or C", then this action takes 5+(R"') onto itself or iS+(C") 
onto itself, as appropriate. Similarly, if T lies in S'L(R") or in ^^(C"), then 
this action takes 7Vf(R") onto itself or A^(C") onto itself, as appropriate. 
The exponential of a linear transformation A on R" or on C" is defined 

by 

oo 

(2) exp^^j:— , 

j=o J- 

and is always an invertible linear transformation on the same space. When 

A is self-adjoint, cxpA is self-adjoint and positive-definite, and indeed the 
exponential defines a one-to-one mapping of iS(R"), iS(C") onto iS_|_(R"), 
5+(C"), respectively. If we write iSo(R"), 5o(C") for the linear subspaces 
of <S(R"), S{C^) consisting of self-adjoint linear transformations with trace 
0, then the exponential is a one-to-one mapping of <So(R"'), <So(C") onto 
A1(R"), A<(C"), respectively. 

Next we look at some of these objects in terms of vector calculus. 

For linear transformations on R" or on C", the differential of the expo- 
nential mapping at the origin is equal to the identity. Explicitly, for each 
linear transformation A on R" or on C", we have that 

(3) f^exp(iA)) =A 



dt 



t=o 



Similarly, for the second differential, if A, B are linear transformations on 
R" or on C", then 



Let us define a Riemannian metric on <S+(R"'), »S+(C") as follows. If P 
is a positive-definite self-adjoint linear transformation on R"^ or on C", then 
tangent vectors to 5+(R'^), 5+(C") at P are given simply by self-adjoint 
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linear transformations A, B on R" or on C*, and we define the inner product 
of these two tangent vectors at P by 

(5) {A,B)p^tT(p-^AP-^B), 

where trC denotes the trace of a hnear transformation C. Using standard 
properties of the trace it is easy to see that this is a positive-definite sym- 
metric bihnear form in A, B, and it depends smoothly on P, so that we get 
a Riemannian metric on 5+(R"^) or on 5+(C"^), as appropriate. 

Suppose that T is an invertible linear mapping on R" or on C", and 
consider the mapping r = T"^ on 5+(R"') or 5+(C"'), as appropriate, defined 

by 

(6) t{P)^T*PT. 

li A, B are self- adjoint linear transformations on R" or on C", which we 
view as tangent vectors to 5+(R"') or iS+(C"') at P, then dTp{A), dTp{B), 
which are the images of A, B as tangent vectors at P under the differential 
of the mapping r, are self-adjoint linear transformations on R" or C" given 

by 

(7) dTp{A) ^T*AT, dTp{B) ^T*BT, 

and they are viewed as tangent vectors to <S+(R") or <S+(C"') at t{P). One 
can check that 

(8) {dTp{A),dTp{B)),^p)^{A,B)p, 

which is to say that the Riemannian metrics on (S+(R"), 5+(C"') are invariant 
under the actions of GLCRJ^), GL(C") that we have defined. 

Of course the usual fiat Riemannian metrics are defined as follows. If 
T is an element of 5(R") or of »S(C"), then again two tangent vectors A, 
B to »S(R") or «S(C") at T are given by self-adjoint linear transformations 
on R" or C", as appropriate, and their inner product in the standard flat 
Riemannian metric is given by the usual inner product, namely 

(9) trAB. 

The flatness of this Riemannian metric is reflected in the fact that it does 
not depend on the point T in the space. 

The metrics that we have defined on 5+(R"), iS+(C") reduce exactly to 
the standard fiat metric at the identity operator /. The differential of the 
exponential function at the zero operator is equal to the identity mapping, 

(10) dexpo(^)=A, 
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as we have noted previously, and thus the standard flat metric on the tangent 
space of iS(R") or i5(C") at agrees with the metric that we have defined 
on 5+(R") or iS+(C"'), respectively, at the identity operator / with respect 
to the correspondence given by the differential of the exponential function. 
We would like to show that in fact this agreement works to another term in 
the Taylor expansion. 

Basically this means that if T is a self-adjoint linear transformation on 
R" or on C", which we think of as being near the origin, and if A, B are two 
self-adjoint linear transformations on the same space, which we think of as 
tangent vectors to >S(R") or 5(0") at T, then the images of A, B under the 
differential of the exponential map at T have inner product with respect to 
the Riemannian metric defined above at the exponential of T which agrees 
to second order with the inner product of A in the standard flat metric on 
5(R") or 5(0"). Explicitly, this means that 

(11) (dexpj,(A),(iexpj,(S))expT 

(12) = tr(exp(-r))(dexp2.(A))(exp(-r))(dexp2^(S)) 
agrees with 

(13) IrAB 

up to terms of order 0(||T|p). This is not difficult to check, using the facts 
that exp(-r) = 7 - T + 0(||r|p), 

(14) dexp^(A) ^A-r\{TA + AT) + 0{\\T\\'), 

and similarly for B. 

As a consequence, if C is any self-adjoint hnear transformation on R" or 
C", then 

(15) exp(iC) 

satisfies the equation for geodesies in the space of positive-definite linear 
transformations at t = 0. Basically this is because the line t C satisfies the 
equation for geodesies in the flat space of self- adjoint linear transformations 
at the origin. Of course we are also using the fact that the standard flat 
metric on the space of self-adjoint linear transformations around agrees 
with the Riemannian metric on the space of positive-definite metrics around 
the identity operator I with respect to the exponential mapping as well as 
they do. 
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In fact, we get that 
(16) 



exp(i C) 



satisfies the equation for geodesies in the space of positive-definite matrices 
for all real numbers t. For a fixed real number to, the statement that this 
curve satisfies the equation for geodesies at to is equivalent to the statement 
that exp((to + 1) C) satisfies the equation for geodesies at t = 0. This is in 
turn equivalent to the statement that T* exp{t C) T satisfies the equation for 
geodesies at t = when T = exp(to C'/2), and this follows from the fact that 
exp(tC) satisfies the equation for geodesies at t = 0, since transformations 
of the form P j^* pj" q^-q isometrics on the spaces of positive- definite 
matrices and therefore preserve geodesies. 

Thus exponentials of straight lines through the origin in the spaces of 
self-adjoint linear transformations give rise to geodesies through the identity 
operator / in the corresponding spaces of positive-definite linear transfor- 
mations. By standard results about uniqueness of initial value problems 
for ordinary differential equations this accounts for all of the geodesies in 
the spaces of positive definite linear transformations which pass through the 
identity operator I. Using mappings of the form P T* PT one can move 
these curves to other places and thereby account for all of the geodesies in 
the spaces of positive-definite linear transformations. 
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